Correlated coherent oscillations in coupled semiconductor charge qubits 
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We study coherent dynamics of two spatially separated electrons in a coupled semiconductor 
double quantum dot (DQD). Coherent oscillations in one DQD are strongly influenced by electronic 
states of the other DQD, or the two electrons simultaneously tunnel in a correlated manner. The 
observed coherent oscillations are interpreted as various two-qubit operations. The results encourage 
searching quantum entanglement in electronic devices. 



Quantum coherence of a single particle and a few 
particles (qubits) has been investigated in various sys- 
tems. Two-qubit unitary operations are key ingredients 
for performing quantum algorithms and correlating mul- 
tiple qubits Typical operations, such as controlled- 
rotation (CROT), which rotates the target qubit state 
conditionally on the control qubit state, and SWAP, 
which swaps quantum states of the two qubits, have been 
demonstrated. However usually one type of operations is 
realized depending on the type of coupling (Ising, Heisen- 
berg, etc.) @, SS, 01- Although other operations can in 
principle be designed in combination with some one-qubit 
operations, it is not a simple task to generate different 
operations under the same coupling. Simple sequences 
for shorter operation time or a smaller number of steps 
have been desired to maintain the coherency of the sys- 
tem 0, 0, IH . In this paper, we report coherent dynam- 
ics for realizing multiple two-qubit operations. Artificial 
and tunable charge qubits fabricated in a semiconductor 
nanostructure are suitable for the purpose, since various 
parameters can be controlled by external gate voltages. 
In the case of Ising (electrostatic dipole) coupling be- 
tween charge qubits fabricated in two double quantum 
dots (DQDs), coherent oscillations of an electron in the 
first qubit can be controlled by the second qubit's state 
(CROT operations) [I flo| | . We find that, even under the 
same coupling, two spatially separated electrons in the 
two DQDs change their locations coherently and collec- 
tively (the correlated coherent oscillations). These coher- 
ent dynamics can be used to design CROT, SWAP and 
other quantum operations, each in a single step. 

We consider a system of two charge qubits, each of 
which possesses an excess electron in an orbital of the 
left or the right dot, \L)j or \R)j, of i-th semiconduc- 
tor DQD (i = 1,2) [111, El El- The electrostatic cou- 



pling between the two qubits in the parallel geometry 
in Fig. 1(a) stabilizes anti-parallel configurations, \LR) 
(= \L)i\R)2) and \RL), rather than parallel ones, \LL) 
and \RR). The coupling is expressed as an Ising- type 
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matrices of the i-th qubit and J (> 0) represents the 
strength of the coupling. The two-qubit system can be 



described by the Hamiltonian 
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where the first term describes the energy offset of the 
single-qubit states and the tunneling coupling energy Aj 
of the z-th qubit. A similar Hamiltonian can be found in 
various physical systems, such as superconducti ng q ubits 
[H and ultra-cold atoms in an optical lattice [3 [HI]. 
Therefore, the following arguments can also be applied 
to those systems. 

Figure 1 (b) illustrates the charge configuration of the 
ground state in the £\ — e 2 plane in the small but finite 
tunneling (Ai = A2 <C J). The boundaries separating 
different charge states represent the resonant conditions, 
where some important quantum operations are expected. 
CROT operations are based on the coherent transitions 
of a target qubit controlled by a control qubit Q. For 
example, a transition between \LR) and \RR) is under- 
stood as rotating the state of the first qubit when the 
second qubit state is |i?) 2 . This transition is expected 
at £\ — J/2 [the right vertical line in Fig. 1(b)] and 
schematically shown in the energy diagram of the upper- 
left inset. Here, \LL) and \RL) should be out of reso- 
nance (separated by J) for the CROT operation. An- 
other CROT operation for the transition between \LL) 
and \RL) appears at E\ = — J/2 [the left vertical line 
in Fig. 1(b) and corresponding energy diagram in the 
lower-left inset]. Similarly, CROT operations when the 
second qubit is the target are expected at e 2 = ±J/2 
(the horizontal lines). 

In contrast, the transition between \LR) and \RL) 
can not be expected as a first-order process, and re- 
quires simultaneous tunneling of two electrons via an- 
other state \RR) (or \LL)) as illustrated by the solid 
(dashed) arrow in the upper-right inset of Fig. 1(b). 
This second-order process takes place when Ai and A 2 
are nonzero. When the first-order transitions are well 
suppressed (||e| ± J/2| ^> Ai, A 2 ), the second-order pro- 
cess dominates at e — E\ = e 2 [on the diagonal line 
in Fig. 1(b)]. In this condition, eigenstates of H 2q 
are approximated to A= {\LR) ± \RL)), \RR) and \LL), 
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FIG. 1: (Color) (a) Colored SEM image of the control device 
(blue and black respectively for unetched and etched surface, 
and gold for metal gates). The circles represent quantum 
dots in DQD1 and DQD2. All measurements were performed 
at electron temperature of 100 mK at magnetic field of 0.6 
T. The energy offsets, Ei and £2, and tunneling coupling, 
Ai and A2, were independently controlled by changing some 
gate voltages simultaneously to compensate for electrostatic 
crosstalk. Tunneling rates are r_Li ~ r_Ri ~ Yl2 ~ T_H2 ~ 1 
GHz for the left (L) and right (R) barrier of the first (1) and 
second (2) DQD. (b) and (c) Charge diagram of the ground 
state at Ai = A2 ~ in (b) and of the first excited state at 
Ai = 13 jueV, A 2 = 25 jueV and J = 25 ^eV in (c). Col- 
ors represent the charge state, LL (cyan), LR (magenta), RL 
(yellow) and RR (white). The resonant conditions are in- 
dicated by solid and dashed lines. Energy diagrams at some 
points (> < 0©) for CROTs, SWAP, and FLIP operations are 
shown in the insets, (d) Electrochemical potential of DQD1 
and DQD2 in the steady state for initialization at large bias 
(Inil and Ini2), in the coherent evolution at zero bias (Evol 
and Evo2), and in the measurement at large bias (Measl and 
Meas2). A schematic of the voltage pulse is shown at the top. 

and the energy gap between (\LR) ± \RL)) is given by 

fl = ±A 2 J/\e 2 - J 2 /4| for equal coupling A = Ai = A 2 . 
This determines the transition frequency between \LR) 
and \RL), and the approximate SWAP operation is ex- 
pected in a single step of the half period. Approxi- 
mate vSWAP operation of the quarter period is useful 
in preparation for correlated states ^= (\LR) ± i\RL)). 

Similarly, another second-order transition between 
\LL) and \RR) [See energy diagram in the lower-right 
inset of Fig. 1(b)] is allowed at e = E\ = —e 2 [on the 
diagonal line in Fig. 1(b)], which can be used to flip 
the total charge polarization of the two-qubit system in 
a single step (We call this process FLIP). The approx- 
imate VFLIP operation is useful in preparation of cor- 
related states A= (\LL) ±i\RR)). In this way, multiple 



two-qubit operations (CROTs, SWAP and FLIP) can be 
performed each in a single step. 

We implemented a two-qubit system by integrating 
two sets of DQDs in a GaAs/AlGaAs heterostructure, 
as shown in the scanning electron micrograph (SEM) in 
Fig. 1(a) [l6|. The two DQDs with individual source 
and drain electrodes are electrically isolated, and thus 
independent currents, I\ and I2, can be measured simul- 
taneously. All qubit parameters can be controlled by 
11 gate voltages. We have already confirmed the elec- 
trostatic coupling between the DQDs from the resonant 
tunneling characteristics, in which the resonant tunnel- 
ing of the first DQD is switched by the charge state of 
the second DQD 0. 

The pulse sequence consists of the following three 
steps, similar to our previous one-qubit experiment [13]. 
First, the system is initialized by setting DQD1 in the 
dissipative single-electron tunneling regime at Vdi = 700 
//V, as illustrated in Inil of Fig. 1(d). An electron is pre- 
pared in the left dot with a high probability (long dwell 
time in the left dot) [13] ■ This initialization works at any 
energy offset e[ (Hereafter, the prime is used for the bi- 
ased situation). On the other hand, DQD2 is kept in the 
Coulomb blockade region with a small bias Vd2 ~ 10 /iV 
(See Ini2), and thus the second qubit is relaxed in the 
steady state. Therefore, the system is initialized in \LL) 
at e' 2 < — J/2 or \LR) at e' 2 > — J/2. Then, the system is 
suddenly brought into the Coulomb blockade condition 
by applying a square (negative) voltage pulse (Vdi = 0) 
with a rise time of about 0.1 ns for a period of t p = 0.08 
- 2 ns, where coherent time evolution is expected (See 
Evol and Evo2 illustrated for the SWAP action). The 
readout of the final state is performed by restoring the 
large bias (Vdi = 700 fiV), where the electron in the right 
dot escapes to the drain and contributes to the current 
(Measl). The above sequence is repeated at 100 MHz 
to obtain measurable current. We obtained the net elec- 
tron numbers flowing per pulse, Npi and Np 2 , for each 
DQD by using lock-in amplifiers [12j. Inelastic tunneling 
during the initialization period (10 ns — t p ) gives a back- 
ground artifact proportional to t p [-0.12 ~ -0.16xi p (ns)], 
which was partially subtracted to highlight the coherent 
oscillations. 

We first demonstrate the CROT operation of DQD1 
using DQD 2 as a control qubit. Figure 2(a) shows the 
coherent oscillations of DQD1 starting from the initial 
state \LR) prepared at e' 2 = +35 /ieV. The overall be- 
havior of the oscillations, including the dependence of 
the period, amplitude, and decoherence time on the de- 
tuning from the resonance (ei = J/2), is consistent with 
our previous study on a single qubit [12j and with the 
two-qubit simulation, from which we obtained Ai = 13 
fieV. When the initial state is prepared to be \LL) at e' 2 
= -60 /iteV, similar oscillations are observed as shown in 
Fig. 2(b) but with a resonance appearing at a signifi- 
cantly different value of £\. Crossover between the two 
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oscillations is depicted in Fig. 2(c), where N p ± measured 
at t p = 0.25 ns (corresponding to 7r-pulse on resonance) 
is plotted in the E\ — £2 plane. The resonant conditions 
discussed in Fig. 1(b) are superimposed in Fig. 2(c) for 
clarity. The vertical patterns of the interference fringes 
are horizontally shifted discontinuously when the initial 
state is altered from \LR) in the upper region to \LL) in 
the lower region. The horizontal shift is a measure of the 
electrostatic coupling J — 25 peV. The large J as com- 
pared to Ai (13 peV) ensures reasonable controllability 
of CROT. 

The above CROT experiments were performed at small 
tunneling coupling of DQD2 (A2 ~ 3 peV estimated from 
an independent measurement). As A 2 is increased while 
keeping Aj constant, the second-order transitions appear 
as additional features (labeled by □) elongated in the 
upper-right direction in Fig. 2(d) for A2 ~ 15 peV and 
in Fig. 2(e) for A2 ~ 25 peV (determined from the fol- 
lowing analysis). The time evolution of Npi (t p ) is also 
investigated at various £\ and £2- Figure 2(f) shows the 



FIG. 2: (Color) (a) and (b) First-order coherent oscillations 
(CROTs) for the initial state \LR) in (a) and for \LL) in (b). 
(c) to (e) Npi in the si — £2 plane obtained at a fixed t v 
= 0.25 ns (-7r-pulse for the CROTs). The energy offsets in 
the initialization period (£1 and e'2) are shown in the right 
and top scales [e[ ~ £1— 50 peV and e' 2 ~ £2 — 20 peV]. 
(f) Correlated coherent oscillations at £1 ~ 25 peV along the 
solid line in (e). (g) Density matrix simulation of (f). (h) 
Typical N p i (t p ) traces; the first-order [ £1 = — J/2 of (b)] 
and correlated coherent oscillations (£2 ~ 10 peV). 



ei dependence of Npi (t p ) at e 2 = 10 ^eV [the dot-dashed 
line in Fig. 2(e)]. Two types of oscillations are resolved: 
one at £\ ~ with the frequency comparable to that of 
the first-order process and the other at E\ ~ 25 peV with 
a lower frequency of 1.3 GHz. Figure 2(h) shows typical 
Npi (t p ) plots of the two cases. We identify the slower 
oscillations to be correlated dynamics of the two qubits 
from the following analysis. 

The resonant conditions (£1,62) were experimentally 
determined at the maximum Np± for sufficiently long 
pulse (t p ~ 1 ns) [for example, the dot-dashed line in 
Fig. 2(f)], and are plotted by crosses (x) in Fig 2(e). 
These conditions can be reproduced just by considering 
the eigenstates of Hiq ■ The resonant conditions were nu- 
merically derived for the minimum energy gap between 
the eigenstates, and are shown by solid and dashed gray 
lines in Fig 2(e) with a fitted parameter A 2 = 25 peV. 
One of these line fits well with the experimental resonant 
conditions. The same resonance lines together with the 
charge state (shown by colors) of the first excited state 
are shown in Fig. 1(c) with the same parameters. As rec- 
ognized by the neighboring charge states, the condition of 
interest (the thick line a — at) represents the resonance 
of \LR) and \RL) (SWAP) for e 2 > J/2, that of \LL) 
and \RR) (FLIP) for £ 2 <C — J/2, and that of superpo- 
sitions of four bases for £2^0. Namely, the straight 
resonant conditions in Fig. 1(b) for Ai = A2 are curved 
as in Fig. 1(c) for unequal coupling (Ai < A 2 ). The ob- 
served coherent dynamics at £2 ~ involves complicated 
superposition of four two-qubit bases, and thus may be 
called correlated coherent oscillations rather than simpler 
SWAP or FLIP operations. 

We also performed density-matrix simulations using 
the standard Lindblad master equation [Tt], EH . In ad- 
dition to the coherent processes described by H 2q , in- 
coherent tunneling transitions to/from the source and 
drain electrodes and spontaneous phonon emission in 
each DQD were included with realistic parameters. To 
reproduce the experimental pulse sequence, the time- 
dependent reduced density matrix p (t) in the Coulomb 
blockade region (at Vp>i = 0) was numerically calculated 
from the initial state p prepared in the transport re- 
gion (at finite Vbi). Expected tunneling electrons (mea- 
surement outcome) in the subsequent transport region 
after the pulse length t p were evaluated. The simulated 
N p i is plotted in Fig. 2(g) by using the same param- 
eter in Fig. 2(f). Although the correlated oscillations 
are significantly degraded by other decoherence mecha- 
nisms (charge noise, finite rise-time of the pulse, etc.), 
the overall oscillation characteristics (fast and slow oscil- 
lations, £i-dependent period around the resonance) are 
well reproduced in the simulations. 

The correlated tunneling can be evidenced by mea- 
suring the signal in the DQD 2. However, no readout 
signal in Np 2 for DQD2 was obtained in the above mea- 
surements, since they were performed in the Coulomb 
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blockade condition of DQD2 [See diagrams in Fig. 1(d)]. 
We confirmed a small but finite negative signal in Np2 
(but positive in iVpi) when DQD2 was made closer to 
the conductive region by adjusting the electrochemical 
potential 62 close to zero (data not shown). This indi- 
cates a fraction of SWAP action (transition from \LR) 
to \RL) produces negative current in DQD2). Since the 
second-order coupling decreases with increasing |ei|, we 
could confirm coherent oscillations in the limited range 
of |ei| < 25 ^cV. Simpler SWAP and FLIP operations 
are expected at larger |ei|. 

We examined such conditions from dc measurement. 
Figure 3 shows (a) I\ at Vjjx = 700 fiV and (b) I 2 at Vjj2 
= 10 /iV. First-order tunneling current in DQD1 and 
DQD2 is seen as broad vertical lines [outside the plot 
range of (a)] and a horizontal line at around £ 2 = ± J/2 
in (b). The second-order cotunncling is simultaneously 
recorded in both currents and appears as very sharp cur- 
rent peaks a, f3, and running in the diagonal direc- 
tions. Peak a is associated with the two-qubit system 
of interest, and (3-/3 involves another orbital state (an 
excited state) neglected in the model. The negative cur- 
rent correlation (positive I\ but negative I2) appearing 
in the upper region (e' 2 > 0) is understood as the cotun- 
neling transition from the initial state \LR) to the final 
state \RL) at e[ ~ e' 2 > J/2 [See arrows in the upper- 
right inset of Fig. 1(b)]. On the other hand, the positive 
correlation (positive I\ and I2) appearing in the lower 
region (e' 2 < 0) of Fig. 3 arises from the transition from 
\LL) to \RR) at e[ ~ -e' 2 > J/2 [See the lower-right 
inset of Fig. 1(b)]. The overall second-order tunneling 
peak constitutes the resonant condition a — al in Fig. 
1 (c) . Although coherent oscillations were confirmed only 
at £2 ~ 0, the appearance of cotunneling peak in the wide 
region in Fig. 3 suggests that coherent SWAP and FLIP 
operations are feasible. 

Finally we'd like to make all data consistent. Data in 



Fig. 3 were taken at the same condition as in Fig. 2(e) 
but in the absence of a pulse (t p — 0). Negative signals, 
where coherent dynamics vanishes, along the diagonal 
lines marked by a, (3, and in Fig. 2(e) coincide with the 
cotunneling currents in Fig. 3, and are attributed to the 
second-order cotunncling occurring in the initialization 
period. This spoils the initialization and measurement 
scheme and thus the coherent dynamics disappears. Such 
experimental details were all reproduced in the density- 
matrix simulation (not shown). 

In summary, we have investigated first- and second- 
order tunneling processes by measuring pulse-induced 
current or current correlations. We expect that multi- 
ple two-qubit operations (CROT, SWAP and FLIP) can 
be induced just by changing parameters £\ and £2- Inter- 
estingly, the two electrons in our device have no spatial 
overlap of wavefunctions, and the expected quantum pro- 
cesses therefore allow us to correlate two electrons (a few 
hundred nanometers apart) without their touching each 
other. This will be an important test of quantum non- 
locality in mesoscopic electron devices and encourages 
the study of charge-based quantum information even un- 
der a constant coupling [7J, ljj|, 2(3]. 
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Y. Tokura, H. Yamaguchi for valuable comments. 



FIG. 3: Correlated cotunneling current in simultaneous mea- 
surement of Ix in (a) and I2 in (b) in the e[ — e' 2 plane mea- 
sured at Vdi = 700 /ieV and Vd2 = 10 /ieV. Solid and dotted 
lines are current traces at e'2 ~ 50 fieV and -10 /ieV, respec- 
tively. Second-order tunneling peaks a, f3, and j3' apear at 
the same marks in Fig. 2(e). 
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